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We study the influence of dissipation on the switching current statistics of moderately damped 
Josephson junctions. Different types of both low- and high- Tc junctions with controlled damping 
are studied. The damping parameter of the junctions is tuned in a wide range by changing temper- 
ature, magnetic field, gate voltage, introducing a ferromagnetic layer or in-situ capacitive shunting. 
A paradoxical collapse of switching current fluctuations occurs with increasing T in all studied 
junctions. The phenomenon critically depends on dissipation in the junction and is explained by 
interplay of two counteracting consequences of thermal fluctuations, which on the one hand assist 
in premature switching into the resistive state and on the other hand help in retrapping back to the 
superconducting state. This is one of the rare examples of anticorrelation between temperature and 
fluctuation amplitude of a physically measurable quantity. 

PACS numbers: 74.40.-|-k, 74.50.-|-r, 74.45.-l-c, 74.72.Hs 



I. INTRODUCTION 

Temperature is a measure of the energy of thermal 
fluctuations. For example, it is well known that noise 
in electronic components or Brownian motion of small 
particles increase with temperature. But does the ampli- 
tude of fluctuations of physical properties always increase 
with T7 Recently a spectacular exception from this rule 
was reported almost simultaneously by three groups of 
researchers P, 0, It was observed that fluctuations 
of the bias current required for switching of a Joseph- 
son junction (JJ) from the superconducting (S) to the 
resistive (i?) state may suddenly collapse (drastically de- 
crease) at elevated T. It was suggested that the para- 
doxical behavior is caused by the fact that temperature 
does not only provide energy for excitation of a system 
from the equilibrium state but also enhance the rate of 
relaxation back to the equilibrium. The latter strongly 
depends on the damping in the system and under cer- 
tain circumstances can reverse the correlation between 
fluctuations and temperature. 

Dissipation plays a crucial role in decay of metastable 
states, which determines dynamics of various physical 
and chemical processes HI. Switching between S 
and R states in JJ's is one of the best studied exam- 
ples of such a decay. The influence of dissipation on 
the switching statistics of JJ's has been intensively stud- 
ied both theoretically 0, @, 0, S and experimentally 
13, ini, E H Q El • The role of dissipation in decoher- 
ence of quantum systems has been widely discussed [H 
and has recently become an important issue for quan- 
tum computing. JJ's are used in several different ways 



in qubit implementations. For example, current biased 



JJ's are employed in phase qubits [17|, where the dissi- 
pation affects relaxation and decoherence in the qubits. 
Furthermore, switching of JJ's is also used for read-out 
of both flux 18] and charge-phase [l^ qubits. 

Here we present an extensive study of dissipation ef- 
fects on the phase dynamics in moderately damped JJ's. 
For this purpose we prepared several types of junc- 
tions with well controlled and tunable damping param- 
eters. In particular, we study low ohmic Nb-Pt-Nb 
Superconductor-Normal metal-Superconductor (SNS) 
junctions, Nb-CuNi-Nb Superconductor-Ferromagnet- 
Superconductor (SFS) junctions with a diluted feromag- 
netic alloys, and Nb-InAs-Nb Superconductor - two di- 
mensional electron gas - Superconductor (S-2DEG-S) 
junctions, as well as Bi2Sr2CaCu208+5 (Bi-2212) high- 
Tc (HTSC) intrinsic JJ's. The influence of dissipation 
on thermal and quantum fluctuations is studied by tun- 
ing damping parameters of the junctions by tempera- 
ture, magnetic field, gate voltage, introducing a ferro- 
magnetic layer or capacitive shunting. The paradoxical 
collapse of switching current fluctuations with increas- 
ing T was observed in all cases. It is shown that the 
collapse temperature critically depends on dissipation in 
the junction. The phenomenon is explained by interplay 
of two conflicting consequences of thermal fluctuations, 
which on one hand assist in premature switching to the 
resistive state and on the other hand help in retrapping 
back to the superconducting state. The conclusions are 
supported by analytical calculations and numerical sim- 
ulations, which are in quantitative agreement with new 
experimental data presented here as well as with those 
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reported in Refs. [E[2,[3. 

The paper is organized as follows. In section II we 
summarize the results of the Resistively and Capacitively 
Shunted Junction (RCSJ) model, required for the analy- 
sis of switching statistics in moderately damped JJ's. In 
Section III we characterize the JJ's studied in this work 
and describe the experimental techniques. In Sec. IV we 
analyze the magnetic field modulation of the switching 
and retrapping currents. This helps to understand the 
origin of hysteresis in the current-voltage characteristics 
(IVC's) and to estimate the damping parameters of our 
JJ's. In Sec. V the main experimental results on the 
switching current statistics are presented. It is shown 
that the collapse of thermal activation (TA) occurs in 
all moderately damped JJ's and that the macroscopic 
quantum tunneling (MQT) phenomena persists even in 
strongly damped SNS-type JJ's. Finally, in Sec. VI we 
discuss the mechanism of the paradoxical collapse of TA 
in moderately damped JJ's and present numeric and an- 
alytic calculations, which clarify the phase dynamics in 
the collapsed state and support our conclusions. 



II. GENERAL RELATIONS 

The damping parameter of JJ's within the RCSJ model 
is characterized by the inverse value of the quality factor 
at zero bias current, Qg. For junctions with sinusoidal 
current-phase relation (CPR) it is given by: 



(1) 



Here uipo = (2e/co/^C)^^^ is the Josephson plasma fre- 
quency at zero bias, R and C are the junction resistance 
and capacitance, respectively, and Ico is the fluctuation- 
free critical current. Determination of Qq is not trivial: 
Ico must be obtained by measurement and extrapolation 
of the switching current statistics [ll|; the effective ca- 



pacitance is not equal to the geometric one because leads 
form a sort of transmission line with finite inductance; for 
Superconductor-Insulator-Superconductor (SIS) tunnel 
junctions the effective resistance is ill-defined. Conflict- 
ing reports exists on what determines the effective damp- 
ing in SIS junctions: the normal resistance the high 



frequency impedance of circuitry the quasiparticle 
resistance [12l ]. or the up-transformed lead impedance 
[2^. Furthermore, the resistance of SIS junctions is 
frequency and bias dependent due to the strong non- 
linearity of the IVC's. 

At finite bias current, the quality factor, Q{I), is given 
by the same Eq.(l) with replaced by the Joseph- 
son plasma frequency at finite bias, which for junctions 
with the sinusoidal CPR is: Wp(/) = ujpo{l — {I / la))"^)^^^ ■ 
This means that the quality factor is bias dependent, 
Q{I) = Qo(l - (J/J^o)^)^/"- For SIS junctions Q{I) can 




FIG. 1: (Color online). Left: the equivalent circuit of the 
RCSJ model. Right: the mechanical analog of the RCSJ 
model: the tilted wash board potential in the energy-phase 
space for / = 0.57co. Arrows indicate three possible particle 
trajectories after thermal escape for different quality factors. 
For the lowest Q the particle get's retrapped in the next po- 
tential well, while for highest Q it will continue to roll down 
the potential, leading to switching of the JJ from the super- 
conducting to the resistive state. 



be further reduced due to frequency dependence of R. 
Therefore, the phase dynamics may change from under- 
damped, Q ^ 1, to overdamped Q ^ 1 as / ^ Ico- 

The electrodynamics of JJ's is equivalent to motion of 
a particle in a tilted wash-board potential formed by su- 
perposition of the periodic Josephson potential and the 
work done by the current source (the tilt), as shown in 
Fig.l. The particle can escape from the potential well as 
a result of MQT or TA process. At low damping the es- 
caped particle will roll down the potential (switch to the 
i?— state). However, if dissipation exceeds the work done 
by the current source it will be retrapped in subsequent 
wells and return to the S*— state, as shown in Fig. 1. 

The TA escape rate from S to R state is determined by 
an Arrhenius law, and for moderate and high damping 
cases it is given by Q: 



F 



TA 



U}p{I) 

at — — exp 

ZTT 



A[/s(/) 



(2) 



Here /SXJs is the potential barrier for 5* — i? switching, 
see Fig. 1, which for the sinusoidal CPR is 



AiJs ^ (4V2/3)£;jo [1 - ///co] 



3/2 



(3) 



where Ej^ = (fi/2e)/co is the Josephson energy. Damp- 
ing enters only into the prefactor of T^a'- 

(1-K1/4Q2)1/2 



at 



1/2Q. 



(4) 



The MQT escape rate can be written in notations of 
Ref. Si as 



-,1/2 



X{Q) exp 



huJr, 



(5) 
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Here 7(T) is the thermal correction with the character- 
istic paraboUc dependence In 7 oc T^; and in the case 
of strong damping xiQ) ^ 27rV3Q^/2[l - Q'^{8ln{2Q) - 
4.428)] and s(Q) ~ 37r[Q + Q-^]. 

The MQT rate can be strongly affected by dissipation 
since the damping dependent factor s{Q) appears un- 
der the exponent in Eq.(l5]). Qualitatively this is due to 
smearing of quantum levels in the wash-board potential 
with decreasing Q. Indeed, the spacing between levels is 
~ hujp, while the level width is h/RC. From Eq.(l) 
it follows that for Q < 1, the width becomes larger than 
the separation between levels, which leads to suppression 
of the MQT. 

In the moderate damping case the crossover between 
MQT and TA occurs at [7,] 



Tmqt{Q) = 



hujp 

27rfc, 



(1 + —Y/^- — 



(6) 



From Eq.® it follows that the MQT-TA crossover tem- 
perature decreases with Q due to general suppression of 
the MQT, as discussed above. 

An analytic expression for the retrapping rate from the 
R to the S state is known only for strongly underdamped 
JJ's Qo > 1 i: 



I — Iro / Ejo 



AC/. 



27r/cBJ 

EjoQl 



R 



ksT 
I -Iro 

IcO 



(7) 
(8) 



Here Iro is the fluctuation-free retrapping current and 
AUr the retrapping (dissipation) barrier. From Eq.® it 
is seen that retrapping, unlike escape, depends strongly 
on damping [l^ > because Qo appears under the exponent 
in Eq.([7|). For underdamped JJ's the Iro is given by 



Iro 



4/, 



cO 



(9) 



This expression is valid for Qo 5^ 3. For smaller Qo 



Iro/Ico - 1.27299-0.31102Qo-0.02965Qo-l-0.01306Qi^. 

(10) 

Eq. (|10p was obtained from interpolation of the numeri- 
cally simulated IVC's within the RCSJ model This 
expression is valid for 0.84 ^ Qo < 3. For Qo < 0.84, the 
IVC's are non-hysteretic, i.e., Iro — IcO- 



Switching statistics 

For a given switching rate Foi from the initial state to 
state 1, the probability of switching within the infinites- 
imal time interval St is given by Poi((5i) = StToi- The 
probability of staying in the state is Vo{St) = 1 — StToi- 



The probability of staying in the state within a fi- 
nite time interval t is then given by the conditional 
probability of not switching during all sub-intervals St: 
Po(t) = lim5t_o(l - StToiY^^* = exp[-roii]. 

Measurements of switching statistics in JJ's are typi- 
cally performed by ramping the bias current at constant 
rate dl/dt. Since the switching rate roi(/) is bias depen- 
dent, the probability of not switching until current / can 
be written as 



^o(/) =exp 



dl/dt 



roi(/)d/ 



Alternatively, it can be written as 



^o(/) = 1- / Poi{I)dI, 



(11) 



(12) 



where Pqi {I) is the probability density of switching from 
to 1 state. By definition Poi{I)SI is the probability of 
switching in the bias interval from I to I+SI. The proba- 
bility density is one of the most important characteristics 
of switching statistics because it directly corresponds to 
experimentally measured switching current histograms. 

Differentiating Eqs. (|ll|12p with respect to / we ob- 
tain: 



Poi 



Po(/) roi(/) 



di 



di/dt 



Poi{I)dI 



(13) 



This equation has a clear physical meaning: the prob- 
ability of measuring the switching event in the current 
interval from I to I -\- SI is the conditional probability 
of switching during the ramping time St — SI / {dl / dt) 
(first term), and the probability that the system has not 
already switched before (second term). The recurrent 
equation p3p is easily solved numerically and couples 
the switching probability density Poi to the switching 
rate Fqi. 

The probability density for switching from S to R state 
follows directly from Eq. p^ . where and 1 are S and R 
states respectively. To obtain the probability density of 
retrapping, Pr, from R to S*— state we should take into 
account that current is now ramped downwards to zero 
and that Pr — a.t I > Ico'- 



PrU) = 



\dl/dt\ 



'IcO 

I - 1^ PR{I)dI 



(14) 



The probability of remaining in the R state and not 
being retraped until current / is 



P«fl(/) = 1 - PR{I)dI, 



(15) 
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TABLE I; Summary of the studied junctions. 



Junction 


Geometry 


/co(M) 




r*(K) 


Qo(T*) 


SFS 


Nb/CuNi 










#la 


70/50nm 


770 


0.24 


- 


- 


#2a 


25/50nm 


178 


0.21 


0.2" 


1.0" 


#2b 


25/50nm 


34 


0.26 


- 


>0.84 


S-2DEG-S 


2DEG wxL 










#la 


lO^m X 500nm 


2.5 


36 


- 


~0.84 


#2b 


40pim X 400nm 


37 


7.5 


0.8 


1.63 


#3a 


lO/xm X 500nm 


2.3 


38.1 




0.88 


#3b 


40/im X 500nm 


7.3 


11.4 


0.1 


0.95 


C-shunted 












#3a 


10/im X 500nm 


2.5 


34 


0.45 


2.4 


#3b 


40fim X 500nm 


7.4 


10.3 


>1 


>3 


SNS 


170 X 88 nm^ 


262 


0.6 




< 0.84 


Bi-2212 


4 X 2.5 ^m^ 


125 


~40 


75 


5.6 



"at H ~ 0.5Oe 



Thermal fluctuations lead to premature switching and 
retrapping. This means that fluctuations tend to de- 
crease the switching current Is with respect to Icq, but 
increase the retrapping current, In, with respect to Irq. 
Therefore, thermal fluctuations help in returning the 
JJ from R to S state. One should also keep in mind 
that since retrapping critically depends on damping and 
damping depends on bias, retrapping can become promi- 
nent at the switching current even for junctions which 
are underdamped at zero bias. 



III. SAMPLES 

Analysis of dissipation effects on phase dynamics re- 
quires junctions with well defined and controlled damping 
parameters. Ideally such junctions should be RCSJ-type 
with bias independent R; and have R much smaller than 
the high frequency impedance of the circuitry Zq ~ lOOO. 
Previous studies of fluctuation phenomena in JJ's were 
performed predominantly on underdamped, Qq ^ 1, SIS 
tunnel junctions, which are not described by the sim- 
ple RCSJ model with constant, Qo, both due to strongly 
non-linear IVC's and considerable shunting by the high 
frequency impedance 0]. This ambiguity does not ex- 
ists for SNS junctions, which are well described by RCSJ 
model with constant R, typically much smaller than Zq. 

Although the quality factor of SNS junctions is ex- 
pected to be constant, verification of this as well as ex- 
act evaluation of Qo is non trivial. Parameters Icq and 
C still need to be independently defined. Furthermore, 
Eq.(l) is valid only for junctions with sinusoidal CPR. 
The CPR in SNS junctions deviates from sinusoidal 
which changes the plasma frequency Wpo and thus affects 
the effective C entering Eq.(l). The effective C is no 
longer equal to the real, explicitly measurable junction 
capacitance. Therefore, quantitative analysis of dissipa- 
tion effects requires a possibility of tuning the quality 
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FIG. 2: (Color online), a) Sketch of the planar Nb-CuNi- 
Nb junction. The junctions were made by cutting Nb/CuNi 
bilayers by FIB. Panels b) and c) show I —V characteristics 
at different T for junctions made from 70/50 nm and 25/50 
nm thick Nb/CuNi bilayers, respectively. The IVC's of both 
junctions are RCSJ-like and exhibit hysteresis at low T. 



factor by at least as many independent parameters as 
the amount of unknown variables in Eq.(l). 

In this work we focus on the analysis of phase dynamics 
in low ohmic SNS-type junctions with moderate damp- 
ing 1 ^ Qo < 10. Emphasis was made on the ability 
to control and tune the damping parameter of the junc- 
tions. Below we describe five different ways used for tun- 
ing and verification of the quality factor of our junctions: 
together with conventional ways of tuning /co by apply- 
ing the magnetic field or changing temperature, we were 
also able to tune /co by applying gate voltage, adding a 
ferromagnetic material into the junction barrier and by 
tuning C by in-situ capacitive shunting. 

The parameters of the studied junctions are summa- 
rized in Table 1. 



Planar SFS (Nb-CuNi-Nb) junctions 

Planar Nb-CuNi-Nb junctions were made by cutting a 
small Nb/Cuo 47Nio.53 bilayer bridge by a Focused Ion 
Beam (FIB) . A sketch of the SFS junction is 
shown in Fig. 2 a). We made junctions from two types 
of Nb/CuNi bilayers with either 70 or 25 nm thick Nb 
layers. The thickness of the CuNi layer was always 50 nm. 
In-plane dimensions of the JJ's, presented here, were the 
same. Details of sample fabrication and characterization 
can be found elsewhere 
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Fig. 2 b) and c) show the IVC's at different T for 
planar SFS junctions made from 70/50 nm and 25/50 
nm Nb/CuNi bilayers, respectively. From Fig. 2 it is 
seen that the IVC's are consistent with the RCSJ model 
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FIG. 3: (Color online). Temperature dependence of the 
switching and retrapping currents for two SFS junctions made 
from the same 25/50 nm Nb/CuNi bilayer at iJ" = 0. It is seen 
that the hysteresis Is > Ir exists at low T for both JJ's. Solid 
and dashed lines represent Ir{T) calculated within the RCSJ 
and self-heating models, respectively. The inset shows the 
size of the hysteresis 1 ~ Ir/Is at the lowest T as a function 
of the dissipated power ai I = Ir. Note that the hysteresis 
is not proportional to the power at retrapping. 



with constant R. Since R Zq ^ lOOfi, shunting by the 
circuitry impedance is negligible also at high frequencies. 

The critical current depends very strongly on the depth 
of the FIB cut. By varying the depth of the cut we 
were able to fabricate JJ's with three orders of magnitude 
difference in Ico [2j|. To the contrary, the resistance of 
junctions remained almost unchanged i? ~ 0.25f2 as seen 
from IVC's in Fig. 2. Since the in-plane geometries of the 



junctions are the same, C and the thermal conductances 
of the junctions are also similar. Therefore, by changing 
the depth of the FIB cut we could vary in a wide range 
the Qo of the junctions by solely affecting Ico and leaving 
all other parameters intact. 

Fig. 3 shows T— dependencies of switching and retrap- 
ping currents for two JJ's on the same chip. Development 
of the hysteresis, Is > Ir, with T is seen. The JJ's had 
good uniformity of the critical current, as follows from 
the clear Fraunhofer modulation of the critical current 
as a function of magnetic field, shown in Fig. 4. 
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FIG. 4: (Color online), a) Magnetic field dependence of the 
switching and retrapping currents for the planar SFS junction 
#2a at T = 37mK. The solid line is the calculated Iro{II) 
within the RCSJ model for Qo{H = 0) = 2.55. b) Is/Ir vs 
H for the same junction. The solid line is the Ico/Iro within 
the RCSJ model for Qo{H = 0) = 2.55. The dashed line 
represents Ico/Ir calculated for the case when the hysteresis 
is caused solely by self-heating (see sec. IV). 



ent magnetic fields along the long side of the J J. Strong 
modulation of the critical current is seen. The IVC's 
are well described by the RCSJ model with constant 

R ~ om < Zq. 

Some JJ's were anodizcd to remove possible shorts 
caused by redeposition of Nb during FIB etching, and to 
further decrease the junction area. The absence of shorts 
and the uniformity of junctions was confirmed by clear 
Fraunhofer modulation of Ic{H), shown in Fig. 6. De- 
tails of the junction fabrication and characterization were 
described in Ref. [26| and will be published elsewhere. 



Nano-sculptured SNS (Nb-Pt-Nb) junctions 



S-2DEG-S junctions 



Nano-scale SNS junctions were made from Nb-Pt-Nb 
trilayers [2^ by three-dimensional FIB sculpturing. The 
thicknesses of bottom and top Nb layers were 225 and 
350 nm, respectively. The thickness of Pt was 30 nm. A 
sketch of the junction is shown in the inset in Fig. 5. 
Nano-fabrication was required both for increasing R and 
decreasing Ico to easily measurable values. 

The main panel in Fig. 5 shows a set of IVC's at 



T = 3.2K for a Nb-Pt-Nb JJ (170 x 



at differ- 



S-2DEG-S junctions with planar geometry were formed 
by two Nb electrodes connected via the 2DEG (InAs) 
27| . Properties {Ico and R) of the JJ's depend on the 
width of Nb-electrodes (either 10, or 40 ^m), the length 
of the 2DEG (either 400 or 500 nm) and the transparency 
of the contact between the Nb and the 2DEG [2§| . A nar- 
row gate electrode was made on top of the 2DEG, forming 
a Josephson field-effect transistor [2^. This provides a 
unique opportunity to tune properties of the JJ's by ap- 
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FIG. 5: (Color online). I — V characteristics at different H 
for a nano-sculptured SNS junctions at T = 3.2K. JVC's 
are RCSJ-like and exhibit hysteresis at low H . Inset shows a 
sketch of the junction. 



plying a gate voltage, Vg 



Details of sample fabrication 



and characterization can be found elsewhere 

Fig. 7 shows a set of JVC's at T = 30mK for different 
Vg. Jt demonstrates that the critical current is increased 
at positive Vg and strongly suppressed at small negative 
Vg. Note that the resistance of the junction starts to 
increase at substantially larger negative Vg < —^V, at 
which the critical current is already strongly suppressed. 
Therefore the JVC's at Vg > —IV are reasonably well 
described by the RCSJ model with a constant R. 

The majority of the switching current measurements 
were performed on the junction (#26) with a wide 
(40/im) and short (400nTO) 2DEG and good transparency 
of the Nb/2DEG interface. For 40/im wide junctions, 
shunting by high frequency impedance was insignificant 
because of the small junction resistance R ~ 7.5 — lOfi. 
Junctions with narrower 2DEG (lO/im) had proportion- 
ally smaller Ico and larger R ~ 30 — 40^1. All JJ's studied 
here had a uniform critical current distribution, as judged 
from periodic Fraunhofer modulations Ic{H), see Fig. 8. 

Fig. 9 shows T— dependencies of Is and Ir for the 
same J J at two magnetic fields, marked by circles in Fig. 
8. The T-dependence of Is and Ir for S-2DEG-S JJ's 
is similar to that of planar SFS JJ's, see Fig. 3. In both 
cases the Ir is T— independent in a wide T— range. 



Bi-2212 intrinsic Josephson junctions 

Intrinsic Josepson junctions (IJJ's) are naturally 
formed between adjacent Cu-0 layers in strongly 
anisotropic HTSC sing le crystals [33|. IJJ's behave as 
SJS-type junctions [3ll Is^l with high Qq, inspite of the 
d-wave symmetry of the order parameter in HTSC. This 
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FIG. 6: (Color online), a) Magnetic field modulation of the 
switching and retrapping currents for the same Nb-Pt-Nb 
junction as in Fig. 5. b) Is/Ir vs H. The lines represent 
calculations within the RCSJ model for Qo{H = 0) = 1.07 
(solid), and in case when the hysteresis is solely caused by 
self-heating (dashed), see sec. IV. It is seen that the capac- 
itive hysteresis within the RCSJ model disappears abruptly 
at certain Ic{H), while the self-heating hysteresis decreases 
gradually with Jc. 



was confirmed by observation of geometric Fiske reso- 
nances 3^, 131 and energy level resolution in the MQT 
experiments on Bi-2212 IJJ 13511 and YBa2Cu307_5 bi- 
epitaxial c-axis (JJJ-like) JJ's [20[. 

IJJ's were made by micro/nano-patterning of small 
mesa structure on top of Bi-2212 single crystals. Here 
we present data for an optimally doped Bi-2212 single 
crystal with Tc — 94.5if . Mesas were cut in two parts by 
FIB to allow true four-probe measurements. Details of 
mesa fabrication can be found elsewhere [36|. Properties 
of our IJJ's were described in detail before [3lJ, i32]. Fig. 
10 shows JVC's at different T for a Bi-2212 mesa. IVC's 
of IJJ's are non- linear and exhibit strong hysteresis below 
Tc- However, at elevated temperatures, 70 < T < 7^, the 
IVC's are almost linear in a small voltage range [3, HH- 
Each mesa contains several stacked IJJ's. Therefore, 
IVC's exhibit a multi-branch structure due to one-by-one 
switching of stacked IJJ's from S to R state. 

Fig. 11 shows the T— dependence of the most probable 
switching current, Ismax, the retrapping current Ir, and 
the fluctuation-free critical current Ico for a sing le IJJ 
from the same Bi-2212 sample as in Fig. 10 37|. The 
Ico was obtained from the analysis of switching statis- 
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FIG. 7; (Color online). The IVC's of a S-2DEG-S junction 
for different gate voltages at T = 30mA". 
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FIG. 8: (Color online). Magnetic field modulation of switch- 
ing and retrapping currents for the S-2DEG-S junction #26 
at T = 30mK. The solid line represents a simulation within 
the RCSJ model for Qo{H = 0) = 2.35. 



tics [36|. Tlie Ico(T) follows the Ambeokar-BaratofF de- 
pendence typical for conventional SIS JJ's [36]. The Ir 
is T— independent at low T and increases with T up to 
~ 85K. Such behavior is also typical for conventional 
SIS JJ's and is attributed to strong T— dependence of 
the low bias quasiparticle resistance, which determines 
the effective dissipation for the retrapping process [s^ . 
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FIG. 9: (Color online). T— dependence of the switching and 
retrapping currents for the S-2DEG-S junction #26 at mag- 
netic fields marked by two of the circles in Fig. 8. The inset 
shows T— dependencies of normalized fluctuation-free critical 
currents Ico{T)/Ico{T — 0) at four H marked in Fig. 8. 



and appears in undcrdampcd JJ's with Qq > 0.84. 

For the case of Bi-2212 IJJ's the hysteresis can be un- 
ambiguously attributed to a large capacitance caused by 
atomic scale separation between electrodes. From mea- 
surements of Fiske step voltages 



the specific capac- 
itance of our IJJ's was estimated as C ~ 68.5/F//zm^. 
Substituting typical parameters of IJJ's [s^l, the criti- 
cal current density Jc{^.2K) ~ IQ^A/cm?] the large bias 
c— axis tunnel resistivity pc — 30 Slcm and the stack- 
ing periodicity s ~ 1.5 rtm, we obtain Qo(4.2if) ~ 20. 
This value will become up to two orders of magnitude 
larger if we use the low bias quasiparticle resistivity at 
T = instead of pc- In any case, IJJ's are strongly 

underdamped, Qo 3> 1, at T <C Tc and also remain un- 
derdamped in practically the whole T— range T < Tc^ as 
seen from Fig. 11. It has been demonstrated that the 
hysteresis Is/Ir in Bi-2212 IJJ's agrees well with the 
calculated Qo using the specific capacitance of IJJ's [sij . 

On the other hand, explanation of the hysteresis in 
SNS-type JJ's is less straightforward. Typically SNS JJ's 
are strongly overdamped, Qo -^1, and the hysteresis is 
caused by either self-heating [seI , non-equilibrium effects 
[iot , or frequency dependent damping [9i] , rather than the 
junctions capacitance. 



Self-heating 



IV. THE ORIGIN OF HYSTERESIS 

Figs. 2-11 show that the IVC's of all the four types of 
JJ's studied here exhibit a hysteresis at low T. According 
to the RCSJ model the hysteresis is related to damping 



It is known that self-heating can cause hysteresis in 
superconducting weak links with negligible C [s^. In 
this case the "retrapping" current simply represents Is 
at the elevated temperature due to power dissipation at 
the resistive branch of the IVC: 
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FIG. 10: (Color online). Four-probe I — V characteristics 
of a Bi-2212 mesa at different temperatures. It is seen that 
the multi-branch structure and the hysteresis persists up to 
~ 3K below Tc. 



Ir^Is{To + AT{Ir))^Is{To) 



dls_ 
dT 



miR))- (16) 



The temperature rise is given by AT{Iii) = PrRui — 
RRthlji, where Rth is the thermal resistance of the junc- 
tion and Pr is the power dissipation at Ir. Thus, Eq. 
(fT6| for Ir^ becomes: 



7fl~/5(ro)[l-a4j, 



(17) 



where a = -{dIs/dT)RRth/Is{To)- The sohition of 
this quadratic equation yields: 



Ir 



^yl + 4aJg - 1 
2als 



(18) 



The dashed lines in Figs. 4 b) and 6 b) show fits 
to experimental Ir{H) in the self- heating model, Eq. 
p8)) . The solid lines in the same figures represent fits 
within the RCSJ model, Eqs. (1,10). /^(i?) for both 
self-heating and RCSJ models were obtained using a 
single fitting parameter {Rth and Qo(H = 0), respec- 
tively), which are unambiguously determined from hys- 
teresis Is/Ir at H = 0. The self-heating and RCSJ 
models provides almost equally good fits to the Ir{H) 
modulation for the SFS junction in Fig. 4. Similarly, 
T— dependences oi Ir, within the two scenarios are practi- 
cally indistinguishable for this sample, as seen from com- 
parison of dashed and sohd fines in Fig. 3. 

For the SNS J J, the RCSJ model provides a better fit to 
the experimental data than the simple self-heating model, 
see Fig. 6b). In experiment and in the RCSJ model the 
hysteresis disappears abruptly at certain H, while self- 
heating is always present and leads to smooth variation 
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FIG. 11: (Color online). Temperature dependence of 
the most probable switching Ismax, retrapping In and 
fluctuation-free critical Ico currents for the same Bi-2212 IJJ. 



of the hysteresis near the minima. On the other hand, a 
better fit can be obtained if we allow some T— dependence 
of Rth- 

A clue to the origin of hysteresis can be obtained from 
comparison of IVC's of JJ's with identical geometry but 
different Is- As described in the previous section, for 
planar SFS JJ's, a minor variation of the FIB-cut depth 
changes Ico by several orders of magnitude without affect- 
ing other characteristics of JJ's (C, R and Rth), as seen 
from Figs. 2 and 3. The inset in Fig. 3 shows the values 
1 - Ir/Is at T ~ 30mK for three SFS JJ's with similar 
geometry on the same chip. Within the self-heating sce- 
nario, Eg. p^ . 1 — Ir/Is would be proportional to the 
power at retrapping, Pr, as shown by the dashed line in 
the inset to Fig. 3. This seems to be true for JJ's with 
large Ico and Pr- However, provided that the hysteresis 
in JJ's with larger Ico is caused by self-heating and Rth of 
all JJ's are the same, there should be no hysteresis due 
to self-heating for the JJ with the smallest Ico in Fig. 
3. This is indicated by the green (lower) dashed line in 
the main panel of Fig. 3, which represents the self- heating 
Ir{T) for the J J with small Is, calculated from Eq. pB]) 
using the parameter a obtained from the fit Ir{T) for the 
J J with larger Is , shown by the black (upper) dashed line 
in Fig. 3. 

Furthermore, experimental Ir(T) for our junctions are 
almost T— independent at low T. This is in contrast to 
the prediction of the sclf-hcating model, Eq. p8)) : Ir cx 

1 /2 

Is , see the black (upper) dashed line in Fig. 3. Note 
that the flatter Ir{T) dependence can not be explained 
by the flatter Ico{T) in comparison to Is{T), because 
within the self-heating model Ir is correlated with Is, 
not Ico- Exactly the same behavior was observed for S- 
2DEG-S junctions, see Fig. 9. 

Thus, we conclude that self- heating does not satisfac- 
tory explain the hysteresis in SFS junctions, although 
it is probably responsible for a considerable part of the 
hysteresis in JJ's with larger Ig. 

For Bi-2212 IJJ's, self-heating was measured directly 
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[41] for the same mesa. The Rth of the mesa ranged froi 
- 70K/mW at T = 4.27^ to - lOK/mW at 80K. Sim 
the dissipated power aX I < Is at the first branch i 
the IVC never exceeded a few ^W, self-heating can I 
excluded as the origin of hysteresis in IJJ's. 



Capacitance 

If the hysteresis were due to finite junction capacitanc 
then magnetic field modulation of //jo(-ff) should be 
unique function of /co(-ff), given by Eqs. (9,10) with fiel 
dependent Qo{H) = Qo{H = 0)[I,o{H) / I.oiH = 0)]i/ 
as follows from Eq.(l). The corresponding IfiQ(H) curv( 
calculated within the RCSJ model are shown by soli 
lines in Figs. 4, 6, and 8. In all cases the agreemei 
with the experimental data is remarkable, considerir 
that there is only one fitting parameter Qo{H = 0) ft 
each curve (indicated in the figures). 

Next we estimate the capacitances that would be r( 
quired for reaching Qo = 1: C[Q(){H = 0) — 1] 35p 
for SFS #2a (Fig.4), - 4pF for SNS (Fig.6) and - 0.2pr 
for S-2DEG-S #2b (Fig.8) JJ's, respectively. Those must 
be compared with the expected geometrical C of the JJ's. 

The overlap capacitance of the SNS junction, Fig.6, 
is small, ~ a few fF, due to small area of the JJ 
O.OlS^m^). The stray capacitance was estimated to be 
of the same order of magnitude. Therefore, the total 
C of this junction is insufficient for observation of the 
hysteresis within the simple RCSJ model. 

The total (stray) capacitance of wide S-2DEG-S JJ's 
with the gate electrode is estimated to be ^ 0.1 — 0.2pF. 
This C can cause a substantial hysteresis in the junc- 
tion #26 with large Ico and may be just sufficient for a 
tiny hysteresis in the other junctions with smaller Ico- 
This conclusion is also supported by observation of un- 
derdamped phase dynamics in those junctions, as will be 
discussed below. 

We argued above that the hysteresis in SFS JJ's can 
not be caused solely by self-heating. But how could the 
huge C ~ S5pF appear in those planar JJ's with the stray 
capacitance in the range of few fF7 To understand this 
we should consider the specific junction geometry, shown 
in Fig. 2 a). The JJ's are made of Nb/CuNi bilayers. 
The CuNi-layer may act as a ground plane for the JJ 
and may create the large overlap capacitance, provided 
there is a certain barrier for electron transport between 
the layers. The transparency of Nb/Cu interfaces, made 
in the same setup, was previously estimated to be ^ 0.4 
[12) S, [3] . The interface transparency between Nb and 
CuNi is expected to be even smaller due to appearance 
of excess interface resistance between normal metals and 
spin-polarized ferromagnets [i^. For typical values of 
the overlap capacitance C ~ 20 — iOfF/^m^, the re- 
quired C ~ 35pF can originate from the bilayer within 
just '--^ 30 — 40/iTO radius from the JJ. An unambiguous 
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FIG. 12; IVC's of two S-2DEG-S junctions on the same chip 
at T — 30mK, H = 0, before and after in-situ C— shunting. 
Inset shows a sketch of the C— shunted junction. 



confirmation of the presence of the large C in our SFS 
JJ's follows also from observation of the underdamped 
phase dynamics, as reported below (Fig. 16). 



In-situ capacitive shunting 

The frustratingly similar behavior of Iji{T,H) within 
self-heating and RCSJ models hinders discrimination be- 
tween self-heating and capacitive origins of the hysteresis. 
To clarify the origin of hysteresis, we fabricated an in-situ 
shunt capacitor, consisting of 300^m wide AI2O3/AI dou- 
ble layer deposited right on top of S-2DEG-S JJ's. The 
sketch of the C— shunted junction is shown in Fig. 12. 
The IVC's for two JJ's from the same chip before and 
after C— shunting are shown in Fig. 12. The length of 
2DEG in those junctions was 500nm, which results in 
considerably smaller Ico than for the JJ#26 with 400nm 
long 2DEG, Fig.8. It is seen that the hysteresis for both 
junctions increased considerably, while R was little af- 
fected by C— shunting [46[ . The total capacitance of the 
shunt was ~ 5pF, much larger than the initial capaci- 
tance of unshunted junctions C ^ 0.1 — 0.2pF. 

It should be emphasized that introduction of the 
C— shunt improve thermal conductance from the junc- 
tions. Indeed, since the sample was placed in vacuum, 
the C— shunt double layer on top of the 2DEG acts as the 
top heat spreading layer 471, |48[ , creating an additional 
heat sinking channel. Thus, the self-heating hysteresis 
must decrease in the C— shunted JJ. Therefore, increase 
of the hysteresis in C— shunted JJ's unambiguously indi- 
cates that the hysteresis is indeed caused by the junction 
capacitance, rather than self-heating. 

To summarize this section, the hysteresis in the IVC's 
of IJJ's are caused solely by the junction capacitance. 
The C— shunted S-2DEG-S JJ's are underdamped with 
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predominantly capacitive hysteresis. For most SNS-type 
junctions the hysteresis is considerably affected by self- 
heating. However, planar SFS and S-2DEG-S junctions 
with short (400 nm) 2DEG are underdamped, Qo ^ ^, 
at low T. Therefore, a substantial part of hysteresis in 
those junctions has a capacitive origin. This conclusion 
is confirmed below by observation of the underdamped 
phase dynamics in those junctions. 



V. COLLAPSE OF THERMAL ACTIVATION 

Measurements of switching and retrapping current 
statistics were made in a carefully shielded dilution re- 
frigerator (sample in vacuum) in a shielded room environ- 
ment. Measurements in the temperature range 1.2-100 K 
were measured in a He'' cryostat (sample in liquid or gas). 
Switching and retrapping currents were measured using 
a standard sample-and-hold technique. All histograms 
were made for 10240 switching events. 



Collapse in Bi-2212 intrinsic Josephson junctions 

A sporadic switching of simultaneously biased stacked 
IJJ's in the Bi-2212 mesa results in rather chaotic switch- 
ing between quasiparticle branches in the JVC, see Fig. 10. 
This makes the analysis of switching statistics quite com- 
plicated. In addition, strong electromagnetic coupling 
of atomic scale stacked IJJ's in the mesa leads to ap- 
pearance of metastable fluxon states 49, |52| and re- 
sults in multiple valued critical current 



a) 

Bi-2212 



4a[5l|. It has 



been reported that switching histograms of IJJ's can be 
very broad and contain multiple maxima [i^, [13, IsH - 
consistent with frustration caused by the presence of 
metastable states in long, strongly coupled stacked JJ's 

In order to avoid the metastable states, we studied 
switching statistics of a single IJ J [H, [s^l . In the stud- 
ied mesa, one of the junctions occasionally had slightly 
smaller Ico (by ~ 20%) than the rest of the IJJ's. Thus, 
we were able to achieve stable switching of this single IJJ, 
while the rest of the IJJ's remained in the S'— state [s^. 

Fig. 13 shows switching histograms for the single IJJ at 
different T. The black solid lines represent simulated his- 
tograms for conventional TA escape, Eqs. (|2ll3|) . at given 
T, and for corresponding junction parameters and experi- 
mental sweeping rates. Detailed analysis of the switching 
histograms can be found elsewhere |36j . The switching 
histograms are perfectly described by the TA theory up 
to T* ~ 75K. However, at higher T the histograms sud- 
denly become narrower. This is clearly seen from Fig. 
13 b) which represents the effective escape temperature 
Tesc obtained from the fit of experimental escape rate 
by the TA expression Eq.([2), with T = T^sc being the 
fitting parameter. For conventional TA, T^sc = T, as 
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FIG. 13: (Color online), a) Switching histograms of a sin- 
gle IJJ at different T. The solid lines represent theoretical 
results for thermal activation at a given T. Note that the ex- 
perimental histograms initially become wider with increasing 
T, but suddenly become narrower and change the shape at 
T* between IIK and 9,1K. b) The effective escape temper- 
ature vs. T extracted from fitting the switching histograms 
to the TA theory. The solid line represents the prediction of 
conventional TA theory, T^sc — T. A sudden collapse of Tesc 
is seen at T* ~ 75K. 



shown by the solid line in Fig. 13 b). A sudden collapse 
of Tesc at T* ~ 75K is clearly seen from Fig. 13 b). We 
also note that the histograms become progressively more 
symmetric at T > T*. 



Collapse in S-2DEG-S junctions 

S-2DEG-S junctions provide a unique opportunity to 
tune Ico , EjQ and Qo by applying the gate voltage Vg , see 
Fig. 7. Fig. 14 a) shows switching current histograms 
at T = 37mK for the S-2DEG-S #36' (identical to #35) 
at different Vg. Panel b) shows that the most probable 
switching current, Ismax, decreases monotonously with 
increasing negative Vg . Panel c) shows the width at half- 
height, A/, versus Ismax- It is seen that initially his- 
tograms are getting wider with increasing negative Vg, 
due to the increase of TA with decreasing Ejq/T. How- 
ever, at < —0.35V a sudden change occurs and AI 
starts to rapidly collapse. 

Fig. 15 a) shows T^sc vs. T for the S-2DEG-S #26 at 
H = 0,2.32,3.05, and 3.66/iT (marked in Fig.8). In all 
cases we can distinguish three T— regions: 

(i) The MQT regime. At low T, Tesc is independent 
of T. Both the TA-MQT crossover temperature, at which 
saturation of TesciT) occurs with decreasing T, and the 
value of Tesc{T — > 0) decrease with if, which leaves no 
doubts that we observe the MQT state 0, S [13, El • Fig. 
15 b) shows the dependence of the Tesc at T = 20mK as a 
function of Icq (H) . Icq for each H was extrapolated from 
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FIG. 14: (Color online), a) Switching histograms of S-2DEG- 
S #3b' at T* = 37mK for different gate voltages Vg. It is seen 
that the height (inversely proportional to the width) of the 
histograms first decreases, but then start to increase with 
increasing negative Vg. b) Dependence of the most probable 
switching current Ismax on the gate voltage, c) The width 
of the histograms A7 vs. Ismax- A sudden collapse of A/ 
occurs at Vg < — 0.35V^. 



the switching histograms. The inset in Fig. 9 demon- 
strates that Ico{T) / Ico{T 0) for diflFerent H collapse 
into one curve, confirming the accuracy of determination 
of IcoiT,H). The solid hne in Fig. 15 b) shows the fit 
to the TA-MQT crossover temperature, Eq.®, taking 
Qo{H = 0) = 2.35, following from the value of hysteresis 
in the IVC's, see Fig. 8. Obviously, the MQT calcu- 
lations are consistent with the previous conclusion that 
this JJ is underdamped at low T and that the hysteresis 
is predominantly caused by the junction capacitance. 

We tried to perform quantum level spectroscopy 
in the MQT state for this J J, but couldn't observe any 
inter level transitions. This implies that the level width 
^ h/{RC) is of the order of level spacing hujp, which in 
turn is caused by a relatively low Qq for this J J. There- 
fore, Fig. 15 confirms the theoretical prediction 0, Q 
that the MQT occurs even in the absence of well defined 
quantum levels in slightly underdamped and overdamped 
JJ's, Qo ^ 1- Both the absolute value and the character- 
istic parabolic shape of Tesc{T) in the MQT state are in 
good agreement with theoretical predictions, as shown in 
Fig. 15 c). The solid lines in Fig. 15 c) show simulated 
Tesc vs T in the MQT state, calculated from Eq.([5|) using 
the experimental conditions for the data in Fig. 15 a). 
It is seen that both the absolute value and the shape of 
simulated Tesc{T) agree with the experimental data. 
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FIG. 15: (Color online), a) Escape temperature vs. T for 
the S-2DEG-S JJ #2& at four magnetic fields (marked by cir- 
cles in Fig. 8). Three T— regions can be distinguished: the 
MQT at low T, the TA at intermediate T, and the collapse 
region at T > T*. b) T^sc in the MQT state at T = 2Qm.K 
as a function of the fiuctuation free IcO, suppressed by apply- 
ing magnetic field. Solid line represents the fit to Eq.© for 
Qo{H = 0) = 2.35. c) Results of MQT simulations, Eq.©, 
for the experimental conditions in panel a). 



(ii) Thermal activation regime. At intermediate T, 
A/ increases in agreement with TA calculations, T^sc — 
T shown by dashed lines in Fig. 15 a). Slightly larger 
inclination of the experimental TesciT) ma,y be due to 
non-sinusoidal CPR in this SNS-type JJ p^. 

(iii) Collapse of thermal activation. At higher T, 
the width of histograms start to rapidly collapse, leading 
to a downturn ofTesciT). The magnetic field dependence 
from Fig. 15 a) reveals that the collapse temperature T* 
decreases quite rapidly with Ico, ie., the collapse occurs 
at lower T in junctions with smaller Qq. 



Collapse in planar SFS junctions 

Fig. 16 shows switching current statistics at different T 
for the planar Nb-CuNi-Nb junction #2a at H ~ 0.5Oe. 
The scales of both axes were kept constant for all his- 
tograms to facilitate direct comparison of the histograms 
at different T. It is seen that the collapse of TA occurs 
at T* ~ 200mK. 
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FIG. 16: (Color online), a) Switching histograms of the pla- 
nar SFS junction #2a at different T. The histograms become 
wider with increasing T up to T = 200mK (lower row) but 
then start to shrink at T ^ 2507717^ (upper row), b) The 
width of histograms A/ vs. T. The AI{T) follows the TA 
behavior up to T* = 200mK but collapses at higher T. 
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FIG. 17: (Color online). Switching histograms of a single IJJ 
a) below and b) above the collapse temperature T* ~ l^K. 
Symbols represent experimental data, dashed lines - TA prob- 
ability density of 5* ^ 7? switching, and dashed-dotted lines - 
probabilities of not being retrapped P„ij, Ea. (|15p . Solid lines 
show the conditional probability density Psr of switching 
without being retrapped. It is seen that close to the collapse 
temperature, the retrapping process becomes significant and 
effectively "cuts-ofF' thermal activation at small bias. Note 
that both the width and the shape of the histograms change 
at T*. Data from Ref. 0]. 



The shape of switching histograms 

Fig. 17 shows switching histograms of the same single 
IJJ as in Fig. 13 just before and after the collapse. It 
demonstrates that not only the width, but also the shape 
of the histogram changes upon the collapse. At T < T* 
the histograms have the characteristic asymmetric shape, 
perfectly consistent with the TA theory [36[ , as shown by 
the black dashed line (coincides with the blue solid line) 
in Fig. 17 a). However, at T > T* histograms become 
narrower and loose the characteristic asymmetric shape, 
as seen from Fig. 17 b). Such a tendency was observed 
for all JJ's, as can be seen from Figs. 13,14,16. As will 
be discussed in sec. VI below, the transformation of the 
histogram shape at T > T* is caused by the interference 
of switching and retrapping processes. 

Effect of C— shunting 

Fig. 18 shows Tesc vs T obtained from switching current 
statistics, for the same S-2DEG-S JJ's as in Fig. 12, before 
and after capacitive shunting. Apparently, C— shunting 
qualitatively changed the phase dynamics of the junc- 
tions, even though it had a minor effect on Jco and R. 



However, C— shunting strongly affected the quality fac- 
tor of the JJ's. From Fig. 18 it is seen that the switching 
statistics of the underdamped C— shunted JJ #3b is well 
described by the TA theory, for which T^sc — T. To the 
contrary, for the unshunted junction, which is just at the 
edge of being overdamped, Qo — 1, the T^sc decreases 
with increasing T almost in the whole T— range. 

Failure of the thermal activation theory in 
moderately damped junctions 

The observed collapse can not be caused by fre- 
quency dependent damping due to shunting by circuitry 
impedance Indeed, we observed the collapse in pla- 
nar SFS junctions with R < 117, for which such shunt- 
ing plays no role. Neither can it be due to T— depen- 
dence of the TA prefactor at in Eq.(|3]). Indeed, damp- 
ing changes only gradually through T* and enters only 
into the (logarithmic) prefactor at of the TA escape rate, 
Eq.(l2]). Gradual variations of at{T) do not cause any 
dramatic variation of the TA escape rate. Moreover, in 
all calculations presented here we did take into account 
the Q{T) dependence of the TA prefactor at, so that Tesc 
must by definition be equivalent to T for the conventional 
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FIG. 18: (Color online). Escape temperature vs T for the 
S-2DEG-S junctions on chip #3 before and after in-situ 
C— shunting. It is seen that the conventional TA behavior, 
Tesc = T, is restored upon increasing Qo after C— shunting. 



TA and the drastic drop in T^sc at T > T*, can not be 
explained within a simple TA scenario. Therefore, the 
observed collapse of switching current fluctuations with 
increasing T represents a dramatic failure of the classi- 
cal TA theory, which was supposed to be valid even for 
overdamped JJ's 0, 01- 



VI. DISCUSSION 

From Figs. 13-16 it is clear that those very different 
JJ's exhibit the same paradoxical collapse of switching 
current fluctuations with increasing T. A very similar 
collapse was observed also in moderately damped SIS 
type Al-AlO^^-Al and Nb-AlO^^-Nb @] junctions and 
squid's. Therefore, the collapse of TA must be a gen- 
eral property of all moderately damped JJ's. 

The dramatic effect of C— shunting on the collapse T* 
clearly shows that damping has a crucial significance for 
the observed phenomenon. From the experimental data 
presented above it is also clear that T* decreases with 
increasing damping and that for overdamped junctions 
T* 0, see the curve for unshunted JJ #36 in Fig. 18. 
The data also shows that the T* is close to the temper- 
ature at which the hysteresis in IVC's vanishes, compare 
Figs. 3 and 16, 9 and 15 a), 11 and 13, which implies 
that retrapping processes may become important in the 
vicinity of the collapse state. 



Influence of retrapping on the switching statistics of 
moderately damped junctions 

The paradoxical collapse of thermal fluctuations and 
the corresponding failure of the conventional TA theory 
in moderately damped JJ's can be explained by the in- 



fluence of retrapping processes on the switching current 
statistics [H, [3, Hi- Indeed, in moderately damped junc- 
tions Ico and I HO are close to each other. As discussed 
in sec. II, increasing T tends to decrease Is and increase 
Ir. Therefore, at sufficiently high T, both switching and 
retrapping events may become possible at the same bias. 
If so, the criterion for measuring the switching event has 
to be reformulated: 

The probability of switching from the S to the R state 
is a conditional probability of switching and not being re- 
trapped back, during the time of experiment. 



(19) 



Here Psr is the probability density of measuring the 
switching event, Ps is the probability density of switch- 
ing, Eq. (fT3|) . and PnR is the probability of not being re- 
trapped Eq. fTS]) . 

Dashed-dotted lines in Fig. 17 a,b) show P„i?,, cal- 
culated for experimental parameters typical for Bi-2212 
IJJ's. The corresponding quality factors are indicated in 
the figures. From Fig. 17 a) it is seen that at T < T* 
the PnR = 1 in the region where Ps > 0, therefore re- 
trapping is insignificant. However, at T > T*, retrapping 
becomes significant at small currents. The resulting con- 
ditional probability density of measuring the switching 
current, Psr, Eq. ()19p . normalized by the total number 
of switching events, is shown by the solid line in Fig. 17 
b). This explains very well both the reduced width and 
the almost symmetric shape of the measured histogram. 



The collapse temperature 

Fig. 19 a) shows the bias dependence of switching 
AUs and retrapping AUr barriers. As was noted in 
Ref. there is always a current Irq < le < IcO at 
which AUs{Ie) = AUR{Ie), so that switching and retrap- 
ping process become equally probable. However, this will 
have an influence on the switching current statistics only 
in the case when the probability density of switching at 
this current is considerable. For the case Qq = 3, shown 
in Fig. 19 a), le/ico — 0-7. However, at low T the major 
part of switching events will occur at Ismax — IcO > le- 
Such situation is seen in Fig. 17 a) for T < T*: the 
probability of not being retrapped FnR — 1 at the most 
probable switching current Ismax and according to the 
criterium, Eq. (|19p . retrapping has no influence on the 
switching statistics. 

Figs. 19 b and c) represent numerical simulations in 
which we intentionally disregarded T— dependencies of 
IcO — Qo and Ejq for simplicity of analysis. Pa- 

rameters were chosen similar to that for the S-2DEG- 
S #26 &t H = 3.66/iT, see Fig. 9. The two consid- 
ered cases correspond to the estimated capacitance of 
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FIG. 19: (Color online), a) Escape and retrapping barrier 
heights as a function of bias current for Qo = 3 and 5. 
b) Numerical simulations of the T-dependencies of the most 
probable switching Ismax and retrapping iRmax currents for 
two values of Qo- Simulations were made for T— independent 
Jco = 20/iyl (dotted line) and parameters typical for the S- 
2DEG-S #2fe junction, c) The simulated width of switching 
A/s and retrapping A/_r histograms disregarding the mutual 
influence of switching and retrapping processes. A/ is the 
resulting width of histograms taking into account switching 
and retrapping. From Figs, b) and c) it is seen that the 
collapse of A J occurs at the condition Ismax — iRmax- 



Since Ismax decreases, while iRmax increases with T, 
switching and retrapping histograms inevitably will over- 
lap at a certain temperature T* . Fig. 19 b and c) clearly 
demonstrates that the collapse of thermal fluctuations of 
the measured switching current A/ occurs at T* and that 
the T* itself strongly depends on Qo- From the simula- 
tions presented in Fig. 19 it is clear that the collapse 
is not caused by a crossover from underdamped to over- 
damped state since Qo was T— independent in this case. 

The collapse temperature can be estimated from the 
system of equations: 



Fta {Is max 

) ~ {dI/dt)/I,o, (20) 

max )=rTA{T*, Ismax). (21) 

Eq. ipn)) states that the J J switches into the R-state 
during the time of the experiment. From Eas. (|2l20p it 
follows that: 



^UsjlSmax) 



In 



O-t'-^pIcO 
2TT{dI/dt) 



Y. 



(22) 



In the measurements presented here Y ~ 24, as seen 
from Fig. 20 a). The parameter Y is weakly (logarithmi- 
cally) dependent on experimental parameters and, there- 
fore, has approximately the same value in different stud- 
ies of switching statistics of JJ's. From Eqs. pi22p we 
obtain the value of the most probable switching current 
I Smax (disregarding retrapping): 



I Smax I IcQ — 1 



(23) 



This dependence is 
19 b) and agrees with 



where Tj = {aV2 E j^) / {iYke] 
shown by the dashed line in Fig 
numerical simulations (squares). 

Similarly, the most probable retrapping current iRmaa 
(disregarding switching) is obtained from Eas. (|7l21l22p : 



the junction #25 (Qo 1.373) and twice the capaci- 
tance (Qo = 1.942), respectively. From Fig. 19 b) it 
is seen that the most probable switching current Ismax 
decreases, while the most probable retrapping current 
iRmax increases with T as a result of thermal fluctua- 
tions. The width of both switching A/g and retrappling 
A//{ histograms continuously increase with T for conven- 
tional TA, as seen from Fig. 19 c). As expected, switch- 
ing histograms are unaffected by the small variation of 
Qo, so that both Ismax and A/g coincide for the two 
values of Qo. To the contrary, retrapping histograms are 
strongly affected by Qo: the lRmax{T) dependence be- 
comes weaker and AIr smaller with increasing Qo. This 
is caused by the increase of the retrappping barrier, AUr 
with Qo as seen from Fig. 19 a). 



iRmax — IrO + IcO\ -p^^ j (24) 

y ^./oVo 

where X = ln[ . f^^-^/— "/f") /-Sffl^] is the log- 

arithm of the ratio of prefactors of TA rettraping and 
switching rates, Eqs. (|2l7p . The factor X is only weakly 
dependent on experimental parameters and in the first 
approximation can be considered constant (or even ne- 
glected). For the case of S-2DEG-S #2fe, AT ~ 3. Red 
dashed hues in Fig. 19 b) represent lRmax{T) calculated 
from Eq.(IM|) with F = 24 and X = 3, which perfectly 
reproduce the simulated lRmax{T) for both Qo values. 

Knowing Ismax{T) and lRmax{T), we can easily obtain 
T* from the condition (cf. Figs. 19 b and c): 
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FIG. 20: (Color online), a) The height of the escape barrier 
at the most probable switching current as a function of T (at 
T = T*): symbols represent experimental data for S-2DEG-S 
#2fe from Fig. 15 a), the solid line corresponds to Ea. (|22|l . b) 
The normalized collapse temperature vs. the quality factor: 
dashed line represents numerical solution of Ea. (|25[) . symbols 
represent experimental data for different JJ's. 



(25) 



A simple analytic estimation of T* can be obtained 
by observing that Ismax{T) is almost linear in a wide 
T— range, as seen from Figs. 19 b) and 11. In this case 
Eg. ipS)) can be approximated as 



ISniax/I*cQ ^ 1 - (iT, 



(26) 



where (3 = 2/ {iTy^Tl'"^ 



To), 1^0 — -^co[l - 

(To/Tj)^/'^/3] and Tq is some characteristic tempera- 
ture ~ Tj. Substituting Eq.JlSl) into Eq. dUD, taking 
a simple approximation for 7^0, Eq.(7), and neglecting 
T— dependence of Ico we obtain a quadratic equation for 
T*, which yields: 



kB{Y + X) 




JO 



2Ql!3E. 
ttQoJ keiY + X) 



- 1 



(27) 

From Eq.jlTl) it follows that kBT*/Ejo strongly de- 
pends on Qo, but is independent of Ejo because T ap- 
pear in Eqs. (|20l24p only in combination T/Ejq (in the 
case of Eq. ([77]) because P ~ I/Ejq). 

The dashed line in Fig. 20 b) represents the numeri- 
cally simulated T* normalized by Ejq and Icq, (left and 
right axes, respectively) as a function of the quality fac- 
tor Qq. It was obtained by numerical solution of Eq. ([25l) . 
without simpHfications used for derivation of Eq. (l27|) . It 
is seen that for overdamped JJ's, Qq < 0.84, T* /Ico 0. 
The T* /Ico continuously grows with increasing Qq > 0.84 



and saturates at ^ 2K/fj,A for strongly underdamped 
JJ's Qo > 1. 

The symbols in Fig. 20 b) represent experimental values 
of T*//co for the JJ's studied in this work. The experi- 
mental data agrees well with the proposed theory (dashed 
line). The simulated values of T*//co are also consistent 
with experimental data for underdamped SIS- type JJ's 
Al-AlO^-Al, T*//co ^ 1 - 3.3K/iiA [1|], and Nb-AlO^- 
Nb, T*/IcQ ~ IK/iiA (with a reservation that those 
measurements were done on SQUID's, which may have 
different activation energies than single JJ's, Eas. l|3l8[ ). 

The T*//co dependence, shown in Fig. 20 b), is al- 
most universal and explains the paradoxical collapse of 
thermal fluctuations, reported in Refs.P, 0, 0|, as well 
as the new data presented here. For example, recovery 
of conventional TA-switching in C— shunted S-2DEG-S 
junctions, see Fig. ??, is caused by the increase of the 
Qo, which according to Fig. 20 b) result in larger T* for 
the same Ico- Similarly, the decrease of Qo due to sup- 
pression of Ico causes the collapse of TA vs. Vg in Fig. 
14 and the decrease of T* with H in Fig. 15. 

From Fig. 20 b) it is seen that in overdamped junctions 
T* — 0, implying that retrapping is crucially affecting 
switching statistics at any T. In this case T^sc{T) and 
AI{T) decrease at all T. We observed such behavior 
for S-2DEG-S junctions with small Ico, and consequently 
small Qo- The tendency of decreasing T* with Ico is 
apparent from Fig 15 a). 

Therefore, observation of the collapse, i.e., a maximum 
of Tesc{T) at T* > 0, is the most unambiguous indication 
of underdamped, Qo > 0.84, state in the studied JJ's. 
The estimation of Qo from the value of T* confirms our 
assessment of junction capacitances, made in sec. IV. 



Phase dynamics in the collapsed state 

The insight into the phase dynamics at T > T* can 
be obtained from Fig. 20 a), in which the dependence of 
AUs{I — Ismax) VS. T IS shown for the case of Fig. 15 
a). The dashed line corresponds to AUs{Ismax)/kBT — 
24.3 ~ Y obtained from simulations presented in Fig. 19 
and demonstrates excellent agreement with the experi- 
ment. The large value of AUs/ksT implies that the JJ 
can escape from the S to the R state only a few times dur- 
ing the time of the experiment. Therefore, the collapse 
is not due to transition into the phase-diffusion state, 
which may also lead to reduction of A J [l3]- Indeed, 
phase diffusion requires repeated escape and retrapping, 
which is only possible for AUs/ksT ^ 1 [^[Hl]. Careful 
measurements of supercurrent branches in the IVC's at 
T > T* did not reveal any dc-voltage down to ~ lOnV 
for S-2DEG-S JJ's and ~ l^iV for IJJ's. Furthermore, 
the IVC's remain hysteretic at T > T*, which is incom- 
patible with the phase diffusion within the RCSJ model 
As can be seen from Fig. 10, the phase diffusion 
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in IJJ's appears only at T > 90K, meaning that all the 
collapse of TA shown in Fig. 13 b) at 75K <T < 85K 
occurs before entering into the phase diffusion state. 

Therefore, in the collapse state the junction makes a 
few very short excursions from the S to the R state during 
the current sweep, before it eventually switches into the 
R state. However, the number of excursions and the 
total excursion time is so small that it does not lead to 
a measurable dc- voltage in the J J. The occurrence of the 
corresponding phase dynamic state, prior to the phase- 
diffusion, has been observed by numerical modelling and 
discussed in Ref. [§]. 



CONCLUSIONS 

We have analyzed the influence of damping on 
the switching current statistics of moderately damped 
Josephson junctions, employing a variety of methods for 
accurate tuning of the damping parameter. A paradox- 
ical collapse of switching current fluctuations with in- 
creasing temperature was observed in various types of 
Josephson junctions [l|,|3,[3, including low-T^ SNS, SFS, 
S-2DEG-S, SIS, and high-Tc intrinsic Josephson junc- 
tions. The unusual phenomenon was explained by an in- 
terplay of two conflicting consequences of thermal fluctu- 
ations, which on the one hand assist in premature switch- 
ing to the i?— state and on the other hand help in retrap- 
ping back into the S*— state. In this case the probability 
of measuring a switching event becomes a conditional 
probability of switching and not being retrapped during 
the time of the experiment. Numerical calculations has 
shown that this model provides a quantitative explana- 
tion of both the value of the collapse temperature T*, 
and the unusual shape of switching histograms in the 
collapsed state. Based on the theoretical analysis, we 
conclude that the collapse represents a very general phe- 
nomenon, which must occur in any underdamped JJ at 
sufficiently high T. 

The collapse of switching current fluctuations in 
Josephson junctions represents an exception from the law 
of increasing of thermal fluctuations with temperature. 
In the studied case, the "failure" of this general law of 
nature, is caused by coexistence of two counteracting pro- 
cesses (switching and retrapping). It should be empha- 
sized that fluctuations for each of the two processes alone 
follow the law and enhance with T in a conventional man- 
ner. It is, however, remarkable and unusual that fluctu- 
ations may cancel each other and lead to reduction of 
thermal fluctuations of a physically measurable quantity. 

Finally, we note that the reduced width of switching 
histograms in the collapsed state of moderately damped 
JJ's may be advantageous for single-shot read-out of su- 
perconducting qubits, which requires accurate discrimi- 
nation of two close current states. 
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